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Abstract 

^ | We propose a simple method to reproduce the 1 Sq pairing proper- 

■ ties of nuclear matter, which are obtained using a sophisticated model, 

by introducing a density-independent cutoff into the relativistic mean 
| field model. This can be applied successfully to the physically relevant 

Q\ ■ density range. 

[ The So pairing gap A in infinite nuclear matter is obtained by solving 

q ' the gap equation, 

s : 

A(p) = -/ v(p,k) , 1 ; =k 2 dk, (1) 

X: 

with v indicating the antisymmetrized matrix elements of the particle-particle 
interaction v pp . One can see from this equation that the physical ingredients 
are the single-particle energies E k and v pp . In sophisticated microscopic 
approaches, the E k are obtained from Brueckner-Hartree-Fock calculations 
with bare N-N interactions, which are fitted to the phase shifts of the N-N 
scatterings in free space. As for v pp , most calculations employ bare inter- 
actions, while some others employ medium-renormalized interactions, such 
as the G matrices. Approaches involving calculations of the former type are 
based on the view that the gap equation itself possesses a mechanism to evade 
strong short-range repulsions, and, accordingly, use of medium-renormalized 
interactions results in a double counting |I|, |^, |3], [|] . The forty-year history of 
non-relativistic studies of the pairing problem || has shown that all the bare 
N-N interactions that are fitted to the phase shifts give almost identical 
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pairing gaps for the 1 So channel. This is because a separable approxima- 
tion [H can be made for the S wave channels in which a virtual ( 1 <S'o) or 
real ( 3 «Si) bound-state pole exists in the T matrices [|7|, and this leads to an 
approximate relation between the pairing gap and the phase shift applicable 
to the low-density region ||. Medium renormalizations are understood to 
cause the gap to become larger because they weaken the short-range repul- 
sion. Irrespective of whether the medium renormalizations are included, the 
particle-hole polarizations should be considered in the next order according 
to a diagrammatic analysis of the gap equation |§, and it is said that they 
act to reduce the gaps [|TU|, |TTJ, |T2]] . 

Another approach to the pairing in nuclear matter is based on the effective 
interactions that are constructed from the beginning to describe finite-density 



systems. An example is the Gogny force [13], which describes the bulk and 



the pairing properties of infinite matter quite well without any cutoffs WA 



and another is represented by various versions of the Skyrme forces, which 



require cutoffs for the description of the pairing ||15|| . From the viewpoint of 



the double counting of the short-range correlation mentioned above, however, 
the adequateness of the use of effective forces in the particle-particle (p-p) 
channel is not evident M . Although this is still an open problem, the Gogny 
force is said to act as a bare force in the p-p channel |L6] . 



Similarly to the studies discussed above, the first relativistic study of the 
pairing in nuclear matter was carried out in 1991 |HJ by adopting a phe- 
nomenological interaction, the relativistic mean field (RMF) model, which 
succeeded in reproducing the bulk properties of the finite-density nuclear 
many-body systems. But the resulting pairing gaps were about three times 
larger than those accepted as standard in the non-relativistic studies. After 
a five-year blank, various attempts to improve this result have begun. These 
attempts can be classified into two groups: The first one employs v pp which 
are consistent with the particle-hole (p-h) channel, i.e. the single-particle 
states [TJ, [B5], ^T[, and the second one employs v pp which are not ex- 



plicitly consistent with the p-h channel |22|, |23[. In addition to these works 
which are based on the single-particle states of the RMF model, there exists 
another [23] which is based on the single-particle states obtained through the 
Dirac-Brueckner-Hartree-Fock (DBHF) calculation [ 25| . We refer to this as 
the third type hereafter. The result that the calculations of the second and 
the third types give almost identical pairing gaps indicates that the pairing 
properties are determined predominantly by the choice of the p-p channel 
interaction, irrespective of the details of the single-particle states. In addi- 
tion, the feature that the obtained gaps are very similar to those given by 
the non-relativistic calculations adopting bare interactions in the p-p chan- 
nel supports this further. As for the first type, a more elaborate calculation, 
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such as one including the N-N polarizations, would be necessary [50]. As 
a complement to this kind of study, however, simpler method suitable for 
realistic applications are also desirable. Examples for which realistic pair- 
ing strengths are indispensable are studies of the crust matter in neutron 
stars and finite open-shell nuclei. In particular, aside from the practical 
successes of the "relativistic" Hartree-Fock-Bogoliubov (HFB) calculations 
implemented by a non-relativistic force [ 2(| , tractable relativistic v pp derived 



from the Lagrangian of the RMF model are needed to keep the concept of 
the HFB calculation. 

The purpose of this paper is to construct a relativistic effective force which 
can be used also in the p-p channel as the Gogny force in the non-relativistic 
calculation. Therefore, first of all, we consider the difference between an 
effective force and a bare force. In Fig. 1(a) the one-boson exchange v pp 
with the coupling constants of the a-u model, which is the simplest version 
of the RMF model, is shown in comparison with the Bonn-B potential |25f . 
which is an example of the relativistic bare N-N interactions. Their shapes 
differ greatly. This is because the former is constructed so as to reproduce 
the saturation property without the short-range correlations, while the lat- 
ter reproduces it in the DBHF calculation which implies them. This leads to 
the characteristic feature of the former that both the small-momentum nega- 
tive off-diagonal matrix elements and the large-momentum positive ones are 
stronger than those for the latter. Both of them enhance the pairing gap, as 
discussed below. 

The momentum integration in Eq. ([I]) should run to infinity when bare 
N-N interactions are adopted. In contrast, there is room to introduce a 
momentum cutoff when we adopt some phenomenological interactions, which 
are meaningful only for small momenta, as the Skyrme force. Evidently the 
assumption of the RMF model that the nucleon is a point particle cannot be 
justified at sufficiently large momenta. Combining this fact with the strong 
cutoff dependence in the momentum region 3-8 fm -1 in Fig. 1(b) suggests 
the possibility to choose a proper cutoff which describes the pairing gap 
quantitatively. Note that the necessity of cutting off the large-momentum 
repulsion in the i> pp derived from the RMF model has also been suggested 



in studies of medium-energy heavy- ion collisions |27j 28] . This is interesting 



in the respect that two different phenomena, which involve large momentum 
transfers, suggest similar cutoffs in the RMF-based p-p interaction. 

In order to describe superfluidity quantitatively, not only the pair wave 
function 

m = 1 - A ^ (2) 

zj{E k -E k¥ y + ^{k) 
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which determines the gap at the Fermi surface 



A(M = ?/ v(k F ,k)(f)(k)k 2 dk, (3) 

47T 







but also its derivative, which determines the coherence length |29 



i 



K { /„•» WVdk j • 1 > 

should be reproduced. The latter quantity measures the spatial size of the 
Cooper pair. In weakly-coupled systems, in which A(k F ) is determined by the 
diagonal matrix element v (k-p, kp) only and £ ^> d (where d is the interparticle 
distance), A(kp) and £ are intimately related to each other. But this does not 
hold for nuclear many-body systems, and the off-diagonal matrix elements 
v(kp,k) play important roles. Therefore here we attempt to find a density- 
independent cutoff A c for the upper bound of the integrals in Eqs. (0) and 
(U) so as to reproduce, in a wide density range, A(k F ) and £ obtained by 
adopting the Bonn-B potential. In other words, we attempt to introduce an 
extra parameter into the o-lq model to fit the pairing properties described 
by a sophisticated model without changing the bulk properties. 

The outline of the numerical calculations is as follows: We start from 



the a-u model with the no-sea approximation, as we confirmed in Ref. |21 
that the Dirac sea effects were negligible. The parameters used are M = 939 
MeV, m a = 550 MeV, m u = 783 MeV, g 2 a = 91.64, and gfc = 136.2 gO 



The calculations were done for symmetric nuclear matter (7 = 4) and pure 
neutron matter (7 = 2). The pairing gap at each momentum is calculated 
by the gap equation (M with A c and the effective mass equation 

M * = M - 4^ r , M * vjk 2 dk . (5) 
ml 2tt 2 Jo ^k 2 + M* 2 k V ; 

Equations ([!]) and @ couple to each other through 
2 1 - f , Ek — Ek F 



v 




(E k -E kF ) 2 + A 2 (k) / 



E k = Vk* + M* 2 + goJ (u°) . (6) 
We search for the value of A c that minimizes 

2 _ 1 ^ 1 ( ^(^F)RMF — A(A;F)Bonn N \ / ^RMF — ^Bonn^ 1 /„n 

X ~2N^\{ A(k F ) Bonn ) + { ^ ) J ' U 
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Here we assume equal weights for A(kp) and £. The single-particle states are 
determined by the a-cu model in both the "RMF" and the "Bonn" cases, as 



in Refs. |22|] and [23]. The summation with respect to kp is taken as kp = 0.2, 
0.3, • • -, 1.2 fm -1 , i.e. N = 11, since we do not anticipate that the present 
method is applicable to the k F ~ CclS6 ; clS discussed later. 

We found that A c = 3.60 fm -1 minimizes \ 2 f° r 7=4. This value 
indicates that not only the small-momentum part, where v(kp, k) < and 
A(k) > 0, but also the large-momentum part, where v(kp, k) > and A(k) < 
0, contribute (see Fig. 1(a)) as pointed out in Refs. ]22| and |23j. The cutoff 
smaller than 2 fm -1 determined in Ref. []B| leads to cutting the "repulsive" 
part completely, and this corresponds to choosing the plateau around 2 fm -1 
in Fig. 1(b), as proposed in Ref. |Tj| in the case of the Skyrme force. The 
present result does not agree with these previous ones. Reference |18 re- 
ports a result different from ours. As discussed in Ref. pTj, there are two 
reasons for this difference. One reason is that they adopted coupling con- 
stants which reproduced the saturation in the Hartree-Fock approximation, 
not the Hartree (so-called MFT) approximation. This leads to larger pairing 
gaps |31 |. The other reason is the difference in the evaluation of the Dirac 
sea effects. 

Figures 2(a) and (b) show how well the a-u model with A c chosen above 
reproduces A(kp) and £, respectively, obtained with the Bonn-B potential 
for symmetric nuclear matter. One can see some deviations between the two 
models both near kp ~ 0.2 fm -1 and kp ~ 1.2 fm" 1 . As for the former, it is 
quite reasonable that the present model based on the mean field picture for 
finite-density systems does not give a good fit. Actually, in such an extremely 
dilute system, the effective-range approximation for the free scattering is 
quite good ||. As for the latter, the deviation results from the fact that the 



superfluid phase in the RMF model, as well as in the Gogny force case p2 



disappears in a nearly A c -independent manner at somewhat larger kp than 
in the Bonn potential case. This at the same time causes the overall peak 
shift of A(kp) to larger kp and makes £ at large kp small. We should note, 
however, that the critical density or kp, where the pairing gap disappears, 
has not been fully discussed yet. The result for pure neutron matter is very 
similar, except that A(kp) is somewhat larger as seen in Fig. 1(b), and the 
superfluid phase survives up to somewhat larger kp, due to larger values of 
M* than in the symmetric matter case. Consequently, the present method 
gives a good fit still for kp ~ 1.2 fm -1 . The density of neutron matter in the 
inner crust of neutron stars corresponds to 0.2 fm -1 <, kp <, 1.3 fm -1 ||. 
Therefore the present simple method covers the greater part of this range. 
In finite nuclei, pairing occurs near the nuclear surface, where the density is 
lower than the saturation point. The present method gives a good description 
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of this region. 

To summarize, we proposed a method to reproduce the 1 Sq pairing prop- 
erties of infinite nuclear matter, obtained using a sophisticated DBHF plus 
a full-range gap equation adopting the Bonn potential, by introducing a mo- 
mentum cutoff into the gap equation with the relativistic mean field model. 
This method was shown to be applicable in a wide and physically relevant 
density range. This points the way to consistent (i.e., using the same inter- 
action in the p-h and the p-p channels) relativistic HFB studies of neutron 
stars and finite nuclei. Finally, we remark that the cutoff we introduced is 
density-independent, since our approach is based on the density-independent 
- except the dependence through M* — RMF interaction which was deter- 
mined at the saturation point. This deserves further investigation. 
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FIGURES 




Figure 1: (a) Matrix element v pp (k, fcp) as functions of the momentum k, 
with the Fermi momentum k F = 0.9 fm" 1 . The solid and dashed curves 
indicate the results obtained by the relativistic mean field model and the 
Bonn-B potential, respectively, (b) Pairing gap at the Fermi surface, k-p = 
0.9 fm _1 , obtained from the relativistic mean field model, as functions of the 
cutoff parameter in the numerical integrations. The solid and dotted curves 
indicate the results for symmetric nuclear matter and pure neutron matter, 
respectively. 




Figure 2: Pairing gap at the Fermi surface (a), and coherence length (b) as 
functions of the Fermi momentum, calculated for symmetric nuclear matter. 
The solid curves and diamonds indicate the results obtained using the rela- 
tivistic mean field model with a momentum cutoff A c = 3.60 fm" 1 and the 
Bonn-B potential, respectively. 
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